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ABSTRACT 

Using flow models based on axisymmetric general relativistic magnetohy- 
drodynamics (GRMHD) simulations, we construct radiative models for Sgr A*. 
Spectral energy distributions that include the effects of thermal synchrotron emis- 
sion and absorption, and Compton scattering, are calculated using a Monte Carlo 
technique. Images are calculated using a ray-tracing scheme. All models are 
scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters 
are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron 
temperature ratio T;/T e . We find that: (1) models with T ; /T e = 1 are inconsistent 
with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly 
increasing function of a*; (3) the X-ray flux is a strongly increasing function of 
i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/T e , but the 
T[/T e = 10 models are generally large and at most marginally consistent with 
the 230 GHz VLBI data; (5) for models with T-jT e = 10 and i = 85 deg the 
event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will 
not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) 
in all models whose SEDs are consistent with observations the event horizon is 
uncloaked at 230 GHz; (7) the models that are most consistent with the observa- 
tions have a* ~ 0.9. We finish with a discussion of the limitations of our model 
and prospects for future improvements. 

Subject headings: accretion, accretion disks — black hole physics — MHD - 
radiative transfer — Galaxy: center - Sgr A* 
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Introduction 



Long term studies of the stellar dynamics in the central parsec of our Galaxy indicate 
that the object in the center of the Milky Way is massive and compact and is therefore 
likely to be a supermassive black hole (we will use Sgr A* to refer to the radio source, the 



putative black hole, and the surrounding accretion flo w). Recent estimates of Sgr A *'s mass 



M = 4.5±0.4xl0 6 M Q and distance/) = 8.4±0.4kpc flGhez et al-lbooal billessen et alkoosl ) 
indicate that it has the largest angular size of any known black hole [GMj {c 2 D) ~ 5.3 /ias). 

Sgr A* is frequently moni tored at all available wavelengths: in radio since its discovery 
in 1974 (iBalick &: Brownll 19741 ). and more recently in submillimeter, near- infrared (NIR), and 
X-rays. It is heavily obscured in the optical and UV (Ay — 30 mag). Sgr A* is a "quiescent" 
galactic nucleus because its bolometric luminosity in units of the Eddington luminosity is 
low, Lb i ~ 10~ 9 L,Edd- The discovery of polarized emission at A = 1.3 mm allowed the use of 
Faraday rotation to place a model de pendent limit on the mass accretion rate 2 x 10 -7 < M < 
2 x 10" 9 M o yr- 1 at r < 20 GM/c 2 faower et allbood . [M7rronegt^alj[2006l). S ubmillimeter 
VLBI of Sgr A* shows structure at very small angular scales (IDoeleman et al.ll2008l ). 



Sgr A*'s spectral energy distribution (SED) can be fi t with semi-analytic 



20051 . lOhsuga et al 



radiatively i nefficient accretion flow (RIAF) models (e.g. iNarayan et al.lll998 
flow models (Yuan et al. 20031). and with time-depe ndent MHD models (e.g. 



quasi-spherical 
), RIAF + out- 



2005 



Goldston et al. 



Moscibrodzka et al 



2006U iHuang et all2007l . Mivoshi et all2008 , 



2007). Other workers have modeled the VLBI 



and submillimeter emiss ion (Broderick & Loeb l2005l.lBroderick fc Loebll2006al. iBroderick fc Loebl 



Broderick et al 



2008 



Yuan et all2009. Dexter et al 



20091 ) assuming a stationary RIAF and computing emission at submillimeter wavelengths 
taking into account general relativistic effects. 

In this work we simultaneously model the spectral energy distribution, including Comp- 
ton scattering, and the VLBI data using a relativistically self-consistent approach. We 
assume that accretion onto Sgr A* proceeds through a geometrically thick, optically thin, 
two-temperature flow that we model using a general relativistic MHD (GRMHD) simula- 
tion. Black hole spin a* is self-consistently accounted for. We also assume that the (likely 
time-dependent, anisotropic, nonthermal) state of the plasma can be described by assigning 
a single temperature T ; to the ions and a possibly different temperature T e to the electrons. 
Conduction is neglected. 

The main goal of this work is to explore how a*, the inclination i, and the ion-to-electron 
temperature ratio T[/T e are constrained by the data. Our paper is organized as follows. In §[2] 
we review broadband observations of Sgr A*. In § [3] we outline our technique for computing 
the evolution of the accretion flow and the emergent radiation. In § H] we present the results 
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of single- and two-temperature SED computations and compare them to the observed SED. 
We summarize and discuss the model limitations in § [5j 



Observations 



Sgr A* has rich observational database in radio ( Serabvn et al. 19971. 



Zhao et al. 2003, An et a. 



2005 



Falcke et al 



Marrone et al. 2006). NIR (Davidson et al. 1992, Herbst et al 



199 



1993UStolow et al 



1996 



elesco et al. 



1996, MentenetalJ 1997, Melia fc Falcke 12001. Hornstein et al 



2002 
2001. 



Genzel et al 



2003 



Belanger et al. 



Eckart et al 



2006, 



2006, 



Aharonianl let al. 20081 ) . 



Schodel et al. 



Baganoff et all hooA JGoldwurm et all I2003L iPorauet et al 



orquet et al. 20081 ) and even 7-rays ( Aharonian et al. 2004 . but 



2007T). X- ravs (Baganoff et al. 



20031. iBelanger et al 



2005 



sec 



In general the emission from Sgr A* in the radio band is rising with the frequency. 
Below v = 10 GHz the s pectral slope a (F u ~ u a ) was found to be a ~ 0.1 (jSerabyn et al. 



1997 



Falcke et al. 



flFalcke et al. 



1998 



199 8]) . Between 10 and 300 GHz the spectral slope changes to a ~ 0.5 



I An et all 120051 ) . 



Marronei (120061 ) reported that the spectral slope becomes flat or declining between 



230 GHz (1.3 mm) and 690 GHz (0.43 mm), consistent with a transition from optically thick 
to optically thin radiation. He estimated a variance-weighted mean value of a = —0.18 from 
four observational epochs (each epoch lasting around 2 hr, and a ch anging from -0.46 t o 
0.08). The quiescent NIR counterpart of Sgr A* has been reported by I Genzel et al.l (120031 ). 
but it is not clear whether the "quiescent" NIR emission a t the position of Sgr A* is back- 
ground noise or a real detection of Sgr A* (IDo et al.l 120091 ) . Thus, the measured quiescent 
emission in the NIR is usually interpreted as an upper limit. The quiescent luminosity at 
the 2-8 keV band measured with the Chandra observatory is Lx = 2.4 x 10 33 ergs _1 , and the 
emission is extended with an intrinsic size of 1.4", consistent with the Bondi radius. The 
mass accretion rate a t the Bondi radius deduced from X-ray observations is ~ lO _6 M yr _1 
( Baganoff et ali 2003 ). Above the stationary emission Sgr A* exhibits intraday variability 
at all observed wavelengths (flares in submillimeter, NIR, and X-rays that often rise simul- 
taneously). 

High frequency VLBI constrains the structure of Sgr A* on angular scales comparable 
to GM/(c 2 D) ~ 5.3/ias. The distribution of intensities on the sky is a convolution of the 



1 Not to be confused with the phenomenological viscosity a of accretion disk theory. In this paper angular 
momentum transport is calculated self-consistently in a GRMHD model. 
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(wavelength depen dent) intrinsic angular structure with ani sotropic interstellar broadening 
proportional to A 2 fjBower et al.ll2006l iDoeleman et al.ll2008l ). Sgr A* has been detected by 
A = 1.3 m m VLBI on baselines b etween Hawaii (JCMT), Arizona (SMTO), and California 
(CARMA) (IDoeleman et al.ll2008l ). This small number of bas elines does not permit imaging 
of the emitting region or the "silhouett e" of the black hole (IBardeenl Il973l . I Chandr asekhar 



19831 . iFalcke et al.ll2000l . lTakahashj|2004j ) . but it does constrain models of the emitting region. 



Using a (two-parameter) symmetric Gaussian brightness distribution model IDoeleman et al. 
f)2008h infer a full width at half maximum FWHM = 37±^//as, or ~ 7GM/(c 2 D). This is very 
small, since the apparent diameter of the black hole is 2^/27GM/ (c 2 D) ~ 55/ias. The FWHM 



for a Gaussian mo del has also been estimated at 7mm (0.268 ± 0.025 mas, or 51GM/( c 2 D), 
Bower et al. l l2004h and at 3.5mm (0.126±0.017 mas, or 28GM/fc 2 £>). IShen et aDl2005h . but 



longer wavelength intrinsic size is more difficult to measure because scatter broadening dom- 
inates the observed image size at A > 1mm. VLBI obs ervations at A = 0.8 m m (345 GHz) 
and 0.65 mm (450 GHz) are expected in the near future (IDoeleman et al.ll2009i ). 



3. Model 

Our model consists of three parts: a physical model of the accreting plasma; a numerical 
realization of the physical model; and a procedure for calculating the emergent radiation from 
the accreting plasma. 

The physical model is a geometrically thick, optically thin, turbulent plasma accreting 
onto a rotating black hole in a statistically steady state. The angular momentum of the 
hole is assumed to be aligned with the angular momentum of the accreting plasma. 
The ions and electrons are assumed to have a thermal distribution fu nction, but with a 



temperature ratio T\/T e that may be different from one (see §2.2.1. in ISharma et al.l 120071 . 
for a discussion of temperature ratios in a collisionless accretion flow model; their work 
suggests that T- x jT c ~ 10 may be a natural value). The equation of state is gas pressure 
P = (7 — l)u (u is the proper internal energy density), with 7 = 13/9, appropriate to a 
plasma with 0, = kTi/{m p c 2 ) <C 1 and O e = kT e /(m e c 2 ) 3> 1 (we will discuss our procedure 
for extracting an electron temperature later). The parameters of the accreting plasma model, 
then, are a* and 7]/T e . 

The numerical realization of the physical model uses the GRMHD code harm, a conserva- 



2 Tilted, or "oblique" accretion flows, require 3D simulations; antialigned flows can be modeled using an 
axisymmetric simulation, but likely provide a worse fit to the data than the low-spin aligned flows considered 
here. 
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tive s hock- capturing scheme with constrained transport to preserve V ■ B = (IGammie et al. 
20031 ). All models in this paper are axisymmetric; we will explore 3D models in a sub- 
sequent publication. Our grid is uniform in modified spherical Kerr-Schild coordinates 
(IGammie et al.l 120031 ). which permit the flow to be followed through the event horizon. 
The coordinates are logarithmic in the Kerr-Schild radius R and nonuniform in Kerr-Schild 
colatitnde 8 (Boyer-Lindquist and Kerr-Schild R and 9 are identical)!] The rotation is 
256 x 256. 



'l-al 



i.e. 



The inner boundary of the computational domain is at Ri n = 0.98(1 + 
just inside the event horizon. The outer boundary is at R out = 40 GM/c 2 ~ 1.8 AU, or 
an angular radius of ~ 210 /xas. Since low frequency emission is believed to arise at larger 
radius, this means that we are unable to model the low frequency (radio and mm) portion 
of the SED. 

We must also supply initial conditions and boundar y conditions. For numerical conve- 
nienc e we adopt the same initial equilibrium torus used by Gammie et al.|[2003 . McKinney fc Gammie 
2004 . and others. This torus has an inner boundary at 6 GM/c 2 and a rest-mass density max- 
imum p max at R max = 12 GM/c 2 Q It is seeded with a weak, purely pol oidal field that f ollows 



the isodensity contours and has minimum (3 = Pj (B 2 / (8tt)) = 100 (see IGammie et al.l (120031 ) 
for details). Small perturbations are added to the int ernal energy. The torus quickly becomes 
turbulent due to the magnetorotational instability (IBalbus &: Hawleylll99ll ). At R < -R max 



the accretion flow soon reaches a nearly (statistically) stationary st ate that is independent 



of the initial conditio ns (except for the magnetic field geometry; see lHawley fc Kroliki 12002 



Beckwith et al.ll2008l ). If our numerical model accurately represents the physical model, this 
inner accretion flow should be similar to the inner portion of a much more extended accretion 
flow. We use outflow boundary conditions at both the inner and outer boundaries, and the 
usual polar boundary conditions at 9 = and 9 = tt. We integrate for 2000GM/c 3 ~ 12 hr, 
or 8 orbital periods at -R max . 



harm (and similar codes) fail if p or u are small in comparison to the kinetic and magnetic 
energy densities, or the density in nearby zones. To prevent this we impose a hard "floor," 

Pmax \ 



so that p > 10~ 4 p max {R/ GM/c 2 )- 3 / 2 and u > lO~ 6 p max c 2 (R/ GM/c 2 )~~ 5/2 . 



To "observe" the numerical model we must specify the observer's distance D and the 
inclination i of the black hole spin to the line of sight. Because the dynamical simulation 



3 The modified Kerr-Schild coordinates xq,xi,X2,X3 are related to spherical Kerr-Schild coordinates 
t, R, 9, by t = xo, R — e Xl , 9 = nx2 + 5(1 — h) sin(27r:E2), and 4> = X3. We set h = 0.3. 

4 For a» = 0.5 we set i? max = 13 GM/c 2 so that scale height H has H / R max similar for all the models. 
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is scale-free but the radiative transfer calculation is not, we need to specify the simulation 
length unit £ = GM/c 2 , time unit T = GM/c 3 , and mass unit Ai (equivalently: the mass 
accretion rate). Since we set the peak density in the GRMHD model to 1 in simulation units, 
the peak density is M. j C? in cgs units. The mass unit is not set by M (which appears in 
the dynamical model only in the combination GM) because the flow mass is <C M, and has 
negligible effect on the gravitational field. M. is therefore a free parameter. 

To calculate the SED we use the relativistic Mo nte Carlo scheme grm onty. A detailed 
description of the algorithm and tests are given in iDolence et al.l (120091 ). The code fully 
accounts for synchrotron emission and absorption, and Compton scattering. It uses a "sta- 
tionary flow" approximation, computing the spectrum through each time slice of simulation 
data as if it were time-independent. This is an approximation because the light crossing 
time is comparable to the dynamical time. It is done because tracking photons through the 
time-dependent simulation data is still too computationally expensive. We will evaluate the 
quality of this approximation once we are able to calculate fully self-consistent spectra. An 
average spectrum is formed by averaging over 50 single slice spectra from each of 4 different 
realizations of the simulation (the realizations differ in the random number seed used to 
generate initial conditions). 

To image the model we use the relativistic ray-tracing code ibothros, which accounts 
onl y for synchrotron e mission and absorption, again using a stationary flow approximation 



I see lNoble et al.l 120071 ). An average image is created using the same averaging procedure as 
for the spectra. 

To sum up, the model parameters (aside from those that describe the initial conditions) 
are M, D, i, a*, T-jT e , and M.. M and D are set by the observations of stellar orbits; we 
briefly explore the consequence of varying them below. A4 we will set for each model by re- 
quiring that time-averaged flux ^(230 GHz) = 3.4 Jy. The remaining three free parameters 
are a*, i, and T-JT C . 

In what follows we explore models with a* = 1 — 2 _n , with n = 1,2,3,4,5 and 6, i.e. 
a, « 0.5, 0.75, 0.88, 0.94, 0.97 and 0.98, T-jT e = 1, 3, and 10, and i = 5, 45 and 85 deg. 

The observational constraints on the model are (aside from the 230 GHz flux den- 
sity) : th^_jubm^meter spectral slope a, the upper limit on the qui e scent NIR flux den 



sitv (IGenzel et al. 



Dodds-Eden et al 



20031 . iMelia fc Falckd 120011 . iHornstein et al.l I2002L iSchodel et all 120071 . 
2009h a nd the upper limit on the quiescent luminosity between 2 and 8 
keV (iBaganoff et al.l 120031 ). Since the measured L x is produced inside 1.4" (« 10 5 GM/c 2 ), 



and our simulation domain covers only the inner 40 GM/c 2 , we exclude models for which 
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Lx is close to the "quiescent emission" Lx- 



4. Results 



We have studied the combinations of model parameters listed in Table [T] (T;/T e =l), 
Tabled (Ti/T e =3), and Table [3] (T-JT C =10). The model that best satisfies the observational 
constraints has a* = 0.94 (model D), T^/T e = 3, and i = 85deg. We will present this "best- 
bet" model in some detail before going on to the full parameter survey (see §4.2p to give the 
reader a physical sense for the models. 



4.1. Best-bet model 



In what follows we will consider time and realization averaged SEDs and images. But 
to get a sense for physical conditions in the accretion flow, consider a single time-slice at 
t = 1680 GM/c 3 . Figure [1] shows the run of n e ,B 2 , and e in the time slice. Evidently 
10 7 cm -3 is a typical equatorial plane density, 30G is a typical field strength, and 6 e = 10 
is a typical electron temperature. Notice that the field changes from a tangled, turbulent 
structure near the midplane to a more organized structure in the "funnel" over the poles of 
the black hole. Temperature generally increases away from the midplane. 

In Figure [2] we show the SED computed from the same timeslice used in Figure [T] (thick 
line) along with the average SED (thin line). The figure also shows a selection of radio, NIR 
and X-ray observational data points (references given in the figure caption). The SED has 
a peak at ~ 690 GHz due to thermal synchrotron emission. Below ~ 100 GHz it fails to fit 
the data because that emission is produced outside the simulation volume. A second peak 
in the far UV is due to the first Compton scattering order, and at > 10 19 Hz (> 40 keV) the 
photons are produced by two or more scatterings. 

Where do the photons originate in each band? Figure [3] maps the points of origin for 
photons in the synchrotron peak (230 to 690 GHz), in the NIR (2 — 10 //m) and in the X-ray 
(2-8 keV). Most of the submillimeter emission originates near the midplane at 4 GM/c 2 < 
r < 6 GM/c 2 . NIR photons are produced in hot regions with high magnetic field intensity 
and high gas temperature, and these are concentrated close to the innermost stable circular 



We use the source brightness profile ( Ba ganoff et al 1 l2003h and estimate that 50% of the X-ray luminosity 
comes from the central pixel of size 0.492". This is still far larger than our computational domain, so we 
require that the X-ray luminosity in our models not exceed Lx = 1.2 x 10 33 ergs _1 . 
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orbit (ISCO; r ISC o ~ 2 GM/c 2 for a* = 0.94) i.e. they come from between 2GM/c 2 < 
r < 6 GM/c 2 . All photons responsible for the formation of the first Compton bump are 
up-scattered between 2 GMj c 2 < r < 8 GMj c 2 but the 2-8 keV emission is produced mainly 
by scatterings in the hottest parts of the disk at 2 GM/c 2 < r < 3 GM/c 2 . Emission can 
be seen around the borders of the funnel in each panel, but this is at a low level and is 
associated with unreliable temperatures assigned by harm to the funnel region. 

The photons that form the submillimeter peak (690 GHz) originate at (r) 4 GM/c 2 , 
where (n e ) »4x 10 6 cm~ 3 (B) « 28 Gauss, and (6 e ) « 17. □ The emissivity vj v {9 = n/2) 
calculated from these averaged values peaks at v ~ 8 x 10 11 Hz, which is quite close to the 
actual peak in vL v (where vL v = ^D 2 i/F v ). The averaged values also yield B 2 /(8tcP) ~ 
0.25; the submillimeter peak plasma is dominated by gas pressure rather than magnetic 
pressure. 

The optically thin synchrotron photons in the energy range between the synchrotron 
peak and first-order Compton scattering bump (6.9 x 10 11 < u < 1.5 X 10 14 Hz) tend to 
arise in current sheets at 2 GM/c 2 < r < 6 GM/c 2 , where 1.5 < n e /(10 6 cm~ 3 ) < 4, (B) = 
30 — 35 Gauss, and e ~ 80 or higher. The corresponding uj u peaks in the mid-IR, at 12 /im. 
These current sheets have higher entropy than the surrounding plasma, consistent with the 
idea that they are heated by (numerical) reconnection. The prominence of the current sheets 
is likely an artifact of axisymmetry. Similar current sheets are seen in axisymmetric shearing 
box models of disks that are not present in three dimensional shearing box models0 

The synchrotron photons Compton scatter at 2 GM/c 2 < r < 8 GM/c 2 where 16 < e < 
70. The X-ray emission at 2-8 keV is formed by scatterings from plasma with (G e ) = 35 
at (r) = 3 GM/c 2 . For this temperature the average energy amplification per scattering is 
1 + 40 e + 160g ~ 1.9 x 10 4 , consistent with the seed photons with energies 2.5 x 10 13 < v < 
10 14 Hz (11 /im > A > 2.9 /mi). This means that many of the seed photons are produced in 
current sheets, and so some uncertainty attaches to the Compton scattered flux. We know 
observationally that Sgr A* produces frequent flares with fluxes larger than those produced 
by our quiescent-source model, so there is a source of seed photons in this energy band, 
albeit a fluctuating one. 

A small fraction of photons are emitted from the funnel wall at large radii (15 — 



6 Angle brackets indicate an average over grid zones weighted with the photon number 'detected' in a 
given frequency band. 

7 Our best guess is that well resolved 3D models will exhibit heating that is more evenly distributed. Direct, 
physical heating of the electrons can tak e place in this strongl y shearing region via plasma instabilities acting 
on anisotropic pressure, as discussed bv lSharma et al.l (|2007l ). 
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40 GM/c 2 ) where the gas temperature is O e ~ 10 3 . This is also likely an artifact of the 
inability of harm and similar codes to track the internal energy of a fluid when the internal 
energy is much smaller than the other energy density scales. Nevertheless, this raises the 
interesting question of what the electron distribution function should be in the funnel. High 
energy electrons might be naturally generated within this tenuous plasma by steepening of 
MHD waves excited by turbulence near the equatorial plane. 

4.2. Parameter survey 

In Figure HI we present averaged spectra for models with different spins (referred to as 
A, B, C, D, E and F; see Tables [fl EJ and [3]), inclination angles i = 85deg, 45deg and 5deg 
in the upper, middle and bottom panels, respectively and temperature ratio T-jT e = 1,3, 
and 10 from left to right. All SEDs are averaged over time and runs as described in § [31 

The tables indicate whether the model is consistent with observations. The model can 
fail in one of four ways: it can produce the wrong submillimeter spectral slope a; it can 
overproduce the quiescent NIR flux; it can overproduce the quiescent X-ray flux; and it can 
be too large at 230 GHz to be consistent with the VLBI data. The last constraint we will 
discuss separately in the next section. It may be useful to recall that M. is adjusted in each 
model so that the 230 GHz flux is 3.4 Jy. 

The model can also fail by cooling too rapidly to be consistent with our neglect of 
cooling in the dynamical model. The Tables list a radiative efficiency r\ = Lbol/Mc 2 , where 
Z/J30L is the bolometric luminosity (integrated over solid angle), and for comparison a thin 
disk efficiency at the same a*, i] ranges between 5.4 x 10 -4 for a* = 0.5, Tj/T e = 10 to 0.18 
for a* = 0.98, T;/T e = 1 (the thin disk efficiency for the latter is 0.25). Only in the a* = 0.98, 
Tj/T e = 1 model is the radiative efficiency sufficiently high that cooling is likely to have a 
significant effect on the GRMHD model. We will consider models with cooling in a future 
publication. 

Very few of the time averaged SEDs based on a single-temperature (T ; /T e = 1) models 
produce the correct a. The exception is edge-on tori {i = 85 deg) around fast spinning black 
holes (model E and F). These models are ruled out, however, because they overproduce NIR 
and X-ray flux. 

For Tj/T e = 3 only the model with a* = 0.94 seen at i = 85 deg agrees with the data. 
This is the best-bet model discussed in § 14. 11 For i = 85 deg, models with spins below 
a* = 0.94 (A, B and C) are ruled out by the inconsistent spectral slope, and models with 
higher spins (E and F), although consistent with the observed a, overproduce the quiescent 



-10- 



NIR and X-ray emission. All models with TJT C = 3 observed at i — 5deg and 45deg are 
ruled out by the inconsistent a. 

For Tj/T e = 10, we find that all models with % = 85deg are ruled out by both a and 
violation of NIR and X-ray limits. For lower inclination angles (i — 5 deg, 45 deg) a few 
models (E and F with % — 5 deg, and A, B, C, and D at % = 45 deg) reproduce the observed 
a. These models are consistent with X-rays and NIR limitations. Models E and F for 
% = 45 deg are ruled out by NIR and X-ray limitations whereas models A, B, C and D for 
i — 5 deg produce a which is too small. 

What is the physical origin of these constraints? 

The dependence on a* arises largely because as a* increases the inner edge of the disk 
- the ISCO - reaches deeper into the gravitational potential of the black hole, where the 
temperature and magnetic field strength are higher. In the disk mid-plane, the temperature 
is a fraction of the virial temperature and scales with radius 6 e oc 1/r. B oc 1/r, while 
the density ~ r, below the pressure maximum. Holding all else constant (which we do not: 
we hold the 230 GHz flux constant) this implies a higher peak frequency for synchrotron 
emission, a constant Thomson depth (in our models the Thomson depth at the ISCO is 
roughly constant, since the path length 1/ ~ risco but the density ~ tisco), and a larger 
energy boost per scattering A ss 160;?, as can be seen in comparing models with different spin 
in Figure 4. The X-ray flux therefore increases with a* because e at the ISCO increases. 

The dependence on Ti/T e is mainly due to synchrotron self- absorption, which is strongest 
at high inclination. For example, because the i = 85 deg, T\/T e = 10 model is optically thick 
at 230 GHz the emission is produced in a synchrotron photosphere well outside risco The 
typical radius of the synchrotron photosphere ranges between 15 GM/c 2 for low spin models 
(a* = 0.5,0.75) and 8 GM/c 2 for high spin models (a* > 0.75). The 230 GHz flux can then 
be produced only with large Ai; as M. increases the optically thin flux in the NIR increases 
due to increasing density and field strength. The scattered spectrum also depends on T\/T e 
since the energy boost per scattering is ~ 160g oc l/(Ti/T c ) 2 . 

The inclination dependence is, interestingly, a relativistic effect. M. is nearly indepen- 
dent of i (it varies by ~ 10%, except for Ti/T e = 10, which due to optical depth effects has 
much larger variation), so models with different inclination are nearly identical. Neverthe- 
less the X-ray flux varies dramatically with i, increasing by almost 2 orders of magnitude 
from i = 5 deg to i = 85 deg. This occurs because Compton scattered photons are beamed 
forward parallel to the orbital motion of the disk gas. The variation of mm flux with i is 
due to self-absorption. The mm flux reflects the temperature and size of the synchrotron 
photosphere. At lower % the visible synchrotron photosphere is hotter than at high %. 
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There is an additional constraint due to Faraday rotation measurements, but this con- 
straint is qualitatively different because we do not directly calculate Faraday rotation in our 
model. Instead we adopt the con straints on M whi c h are inferred, yia a separate model, 
from the Faraday rotation data (IBower et al.l 120051 . iMarrone et al.l 120061 ). M increases, 
in a nonlinear way, with increasing T\/T e . For T;/T e = 1,M > 6 x 10~ 10 M Q yr _1 . For 
Ti/T e = 10, M < 4 x 10~ 7 M© yr _1 . All these values are consistent with the Faraday rotation 
constraints, although the highest M, T- Y /T e = 10, models are only marginally consistent. 

There are a few other general trends worth mentioning. In all models the average 
optical depth drops below 1 at 0.4 to 1.3 mm. For i > 45deg and high BH spins a* > 0.75 
the emission in NIR band (2 /im) is formed by the direct synchrotron emission while the 
2 — 8keV emission results from a first-order scattering. For low BH spins a* < 0.75 the 
emission in NIR is due to first-order Compton scatterings and the X-ray is second-order 
scattering. For T\/T e = 10 and % = 5deg independently of the BH spin the NIR emission is 
formed by a first-order Compton scatterings and X-rays- by second-order scatterings. 



4.3. Images and the size of the emitting region 

We compute the 230 GHz intensity maps of our models using a ray tracing code 



(INoble et al.ll2007l ) and we average them in the same manner as the spectra. To estimate the 
size of the emitting region we calculate the eigenvalues of the matrix formed by taking the 
second angular moments of the image on the sky (i.e. the length of the "principal axes"). 
The eigenvalues along the major (01) and minor (02) axis are given in Table HI In Figure [5] 
we show averaged 230 GHz images for models with SEDs that are consistent with the data. 

The source size depends on the model parameters. For i > 45 deg, we find a critical mass 
accretion rate M ~ 10~ 8 M Q yr^ 1 (the exact value depends on a* and i) below which the size of 
the emitting region decreases monotonically with increasing T ; /T e . Above M « 10~ 8 M Q yr^ 1 , 
the size of the emitting region increases with increasing T±/T e . The increasing trend can be 
explained by the appearance of the synchrotron photosphere at 230 GHz for larger mass 
accretion rates at high inclination. For 2 = 85 deg and Ti/T e = 10, at 230 GHz the black hole 
horizon is cloaked by the photosphere and cannot be observed by VLBI (notice that this 
model is ruled out for other reasons). For a constant T-jT e and i the size always decreases 
with increasing a*, because the emissivity of the central regions increases with a*. 

The size of the emitting region for our best-bet model is consistent with the observed 
FWHM = 37/ms (inferred from VLBI data using a two-parameter Gaussian model). For 
Ti/T e = 10 the sizes of the images are inconsistent with the VLBI measurement, except model 



-12- 



D, which is only marginally consistent. Notice that this moment-based analysis is crude. It 
would be better to "observe" the model with the same baselines used in gathering the VLBI 
data (this would add a new parameter, the position angle). Our analysis is particularly ill 
suited to low i models that are ring-like and therefore poorly fit by a Gaussian model. 



4.4. Varying distance and mass 

In our discussion we have fixed the mass and distance of Sgr A*, but these are uncertain 
to ~ 5%. How would changing these parameters change our results? 

First, consider how M. depends on M and D. Near the submillimeter peak, F v oc 
n e BC 3 D- 2 ~ M z/2 C- l ' 2 T- l D- 2 ~ {M/Mf/ 2 D~ 2 (since C ~ M and T ~ M) if the model 
is optically thin (although there are usually optically thick lines of sight through the model 
even if the mean optical depth at the submillimeter peak is < 1). We therefore expect that 
M ~ D*/ 3 M. 

Consider varying M and D in our best-bet model (model C; T;/T e = 3 and i = 85 deg). 
We find that A.A4/A4 ~ 15% when changing the distance from 8.0kpc to 8.8kpc if we fix 
M = 4.5 x 10 6 M Q . In particular, M = 17.5 x 10 18 , 18.9 x 10 18 , 20.6 x 10 18 for D = 8.0, 8.4 
and 8.8 kpc, respectively. For D = 8.4 kpc, M = 18.2 x 10 18 , 18.9 x 10 18 and 20.6 x 10 18 
for M = 4.1,4.5 and 4.9 x 1O 6 M respectively, which gives AA4/Ai pa 10% when changing 
black hole mass from 4.1 — 4.8 x 1O 6 M . This is crudely consistent with our expectations 
based on an optically thin source. 

Finally we change the mass and distance simultaneously according to the observational 
relation M J D" L8 =constant Jdiez et al.ll2008h . For D = 8.0 kpc and M = 4.1 x 1O 6 M , 



M = 16.6 x 10 18 whereas for D = 8.8 kpc and M = 4.8 x 1O 6 M , M = 20 .4 x 10 18 . We find 
that at D = 8.0 kpc, M = 4.1 x 1O 6 M , a = -0.32, log 10 L x = 33.0 ; at D = 8.4 kpc, M = 
4.5 x 1O 6 M a = -0.44, log 10 L x = 32.9; at D = 8.8 kpc, M = 4.8 x 1O 6 M , a = -0.47, 
log 10 Lx = 32.7. The spectral slope and X-ray luminosity therefore vary < 25%. 

Our best-bet model remains consistent with the data, then, if we vary with M and D 
within the range permitted by observation. Models with T;/T e = 10 and i = 45 deg with 
BH spin a* = 0.97 and 0.98 become acceptable if M and D are lowered, but only the model 
with a* = 0.98 would be (marginally) consistent with VLBI measurements of the Sgr A*size. 
In sum, D and M are tightly constrained; varying them within the narrow range of values 
permitted by observations does not change the main conclusions of this work. 



13 



5. Summary 



Under the assumption that the accretion flow at the galactic center is optically thin, 
geometrically thick, and lightly magnetized, we have presented constraints on a*, T-jT e , and 
i for Sgr A*. We find that models with T ; /T e = 3 and 10 describe the sub- mm spectral 
observations (a) better than models with T-jT e = 1. We find that the model with TjT e = 3, 
black hole spin a* pa 0.94 and the close to edge-on inclination angles is consistent with the 
broadband SED observational data and the size of Sgr A* measured by VLBI. In this case 
the silhouette of the black hole is difficult to observe because Doppler boosting of the disk 
emission places almost all the emission on one side of the black hole. 

If, on the other hand, the electrons are heated relatively inefficiently {T\/T c = 10) then 
models with a* = 0.97,0.98 observed at i — 5deg, or a* = 0.5,0.75,0.88,0.94 observed at 
i = 45deg are consistent with the observed SED. The sizes of the emitting regions in these 
models, however, seem to be inconsistent with the VLBI measurements, except again at 
a* = 0.94. 



Our best-bet estimate of the b lack hole spin 



= 0.94) disagrees with iBroderick et al. 
(|2008|) (following lYuan et all 120091 ) who found a* = +0A and i = 90deg_ 

-sodeg (2c errors) 

based on a careful analysis of images of RIAF models. The discrepancy may be a consequence 
of different emissivity (ours is based on Leung et al. 2009, in preparation), and different 
underlying models for the run of temperature, density profile, magnetic fie ld strength, and 



geome try of the flow. We also do not include non-thermal emission as in IBroderick et al. 
( 120081 ). The results for T;/T e = 10 and % = 45 deg at low spin values agree with the previous 



study, but according to our moment analysis these models are inconsistent with the VLBI 
data. An analysis of images at v — 230 GHz that folds the models through the VLBI 
observation process is needed to definitely exclude models based on the VLBI data. 



Our a* and i constraints are different from those presented in iNoble et al.l ( 120071 ). be- 
cause here we allow T\jT c ^ 1. We do not find a good fit to the observational data for 
single-temperature models which were studied in the earlier work. We confirm the trend 
that the bolometric luminosity increases with the increasing black hole spin. 



The models studied here differ in many respects from those considered earlier by lMoscibrodzka et al 



( 120071 ). The earlier models were based on low angular momentum, nonrelativistic hydrody- 
namic models for the accretion flow that extended over a wide range in radii. The models 
described here are fully relativistic MHD models that extend over a limited range in radius 
and use fully relativistic radiative transfer. 



There are still significant uncertainties in our models. These uncertainties fall into four 
categories: the unimportant; those which may be important and be easily eliminated with a 
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small additional effort; those which may be important and require a major effort, but are in 
principle straightforward to eliminate; and those which are serious and require new physical 
understanding. 

In the interests of full disclosure, our unimportant approximations are: (1) we use a 7 
law equation of state rather than a Synge-type equation of state that would more accurately 
represent our two-temperature relativistic gas. Shiokawa et al. 2009 (in preparation) show 
that the associated changes in spectra are small; (2) for T-JT C > 10 our cyclo-synchrotron 
emissivity and absorptivity are imperfect because most of the emission comes from electrons 
with e ~ 1, where our approximate expression begins to break down (Leung et al. 2009, 
in preparation) 0; (3) we neglect bremsstrahlung, which is expected to be important only 
far from the horizon; (4) we neglect double Compton scattering. The cross section for the 
double Compton process is ~ e 2 /hc (fine-structure constant) smaller than single Compton 
scattering and can be neglected here (also hu <C m e c 2 ); (5) we neglect induced Compton 
scattering. Induced Compton scattering is important for (kBTb/m e c 2 )rT > 1, where TJ, is 
the brightness temperature. In Sgr A* pa 10 10 K and tt <C 1 so it is indeed negligible. 

The significant approximations that could be fixed with some additional effort include: 
(1) our neglect of cooling. It is straightforward to run our GRMHD models with cooling, 
but then they are no longer scale-free. With cooling turned on we would need to fix hA by 
evolving the GRMHD model at a trial Ai, calculating the spectrum, and repeating until 
1*^(230 GHz) = 3.4 Jy. (2) axisymmetry. Three dimensional (3D) models are available but 
far more expensive to evolve. Use of a 3D models would permit us to evolve models with (3) 
a wider range of radii, so that millimeter emission could be included and even the sub-mm 
emission could be more accurately modeled. We have not run axisymmetric models with 
radially extended accretion flows because they tend to develop pathologies (strong, radially 
extended magnetic filaments). (4) our neglect of nonthermal electrons. These could be 
readily included using a phenomenological prescription f or the shape and amplitude of the 



nonthermal porti on of the electron distribution function (jOzel et al.ll2000l . lYuan et al.l 12003 



Chan et al. 20091 ). (5) our use of the steady-state approximation in calculating SEDs and 



images. This would require time- dependent radiative transfer, which is straightforward in 
principle but computationally expensive. 

Other approximations can be fixed only with significant additional effort: (1) our neglect 



8 For models with Ti/T e = 1 and 3 the error in SED associated with our approximate emissivity function 
is less than 1% at all frequencies. For T-JT C = 10 the errors are less than 10% at 230 GHz because the 
emissivity- weighted mean temperature is lower and self-absorption is important. The errors are taken from 
comparison of Leung et al. 2009 emissivity formula with directly integrated cyclo-synchrotron harmonics for 
lowest values of Q e found in our simulations. 
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of pair production. This would require a model for the radiation field near the pair production 
threshold. Our preliminary estimates suggest that in many of our models pair production 
is substantial. One advantage of incorporating pair production is that it might permit us 
to eliminate our numerical floor and therefore more accurately evolve the low density funnel 
region; (2) our treatment of thermal energy in the funnel, harm tends to produce high 
temperatures in the tenuous funnel plasma, some of which are clearly numerical artifacts 
caused by application of the density floor and other, more subtle, numerical issues associated 
with the small ratio of thermal energy density to other energy densities; (3) our simplistic, 
two-temperature thermal model for the plasma. This includes our neglect of conduction and 
anisotropy of the plasma. 

New physical understanding would be required to predictively model (1) nonthermal 
parts of the distribution function and (2) the initial magnetic field configuration. Nonther- 
mal particles can of course be included in a phenomenological prescription, but the particle 
injection and acceleration processes are still not fully understand. As we have already men- 
tioned, prior work shows that the GRMHD models depend nontrivially on the initial field 
configuration. We have adopted a simple, numerically appealing initial configuration, but 
the long-term evolution of the large-scale field is ill understood. 

Finally, notice that there are observational constraints from polarization data and from 
light curves (statistically, the one and two-point statistics of the light curves at each fre- 
quency, and the cross correlations between different frequencies). Treating the polarization 
data requires accurate emissivities and absorptivities, as well as models that extend well 
past the radius where O e = 1, which is where most of the intrinsic Faraday rotation occurs. 
The light curves require full, time-dependent radiative transfer, since the dynamical time is 
comparable to the light crossing time. 

This work was supported by the National Science Foundation under grants AST 00- 
93091, PHY 02-05155, and AST 07-09246, through TeraGrid resources provided by NCSA 
and TACC, and by a Richard and Margaret Romano Professorial scholarship, a Sony faculty 
fellowship, and a University Scholar appointment to CFG. The authors are grateful to Stu 
Shapiro, Fred Lamb, Dan Marrone, Shep Doeleman, and Vincent Fish for discussion and 
comments. 
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Table 1. Summary of MHD and SED for models with T-jT e = 1 



run 


a* 


t 

[degj 


[10 M0 yr ] 


a 


logio Lx 
[ergs -1 ] 


V 


Vtd 


consistent 
w/ obs.? 






5 


3.97 


-2.04 


30.9 


8.9 x 10" 3 




NO 


A 


0.5 


45 


3.55 


-1.78 


30.7 


8.0 x 10~ 3 


0.0821 


NO 






85 


3.81 


-1.21 


31.3 


8.6 x 10" 3 




NO 






5 


1.6 


-1.94 


31.1 


2.4 x 10~ 2 




NO 


B 


0.75 


45 


1.39 


-1.69 


31.0 


2.1 x 10~ 2 


0.112 


NO 






85 


1.5 


-1.12 


31.6 


2.2 x 10~ 2 




NO 






5 


0.94 


-1.92 


31.3 


4.8 x 10~ 2 




NO 


C 


0.875 


45 


0.87 


-1.65 


31.3 


4.3 x 10~ 2 


0.145 


NO 






85 


0.91 


-1.07 


32.3 


4.5 x 10- 2 




NO 






5 


0.86 


-1.67 


31.9 


7.1 x 10- 2 




NO 


D 


0.94 


45 


0.78 


-1.41 


32.1 


6.4 x 10~ 2 


0.179 


NO 






85 


0.80 


-0.87 


33.1 


6.6 x 10~ 2 




NO 






5 


0.85 


-1.51 


32.3 


0.11 




NO 


E 


0.97 


45 


0.78 


-1.25 


32.7 


0.1 


0.213 


NO 






85 


0.80 


-0.68 


33.8 


0.1 




NO 






5 


0.64 


-1.51 


32.9 


0.19 




NO 


F 


0.98 5 


45 


0.58 


-1.25 


33.1 


0.17 


0.245 


NO 






85 


0.59 


-0.69 


34.2 


0.18 




NO 



Note. — The columns from left to right are: run ID, dimensionless spin of the black hole, inclination 
angle of the observer with respect to the black hole spin axis, averaged rest mass accretion rate, a 
spectral slope between 230-690 GHz(F ~ and luminosity in the X-rays (at v ~ 10 18 Hz), the 

radiative efficiency rj = Lbol/Mc 2 , the thin disk efficiency for the same a*, and whether the model is 
consistent with the data. 
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Table 2. Summary of MHD and SED for models with T\/T e = 3 



run 


a. 


i 

[deg] 


(M> 
[lO^Meyr- 1 ] 


a 


log 10 L x 
[ergs -1 ] 


<7 


VTD 


consistent 
w/ obs.? 






5 


10.0 


-2.09 


30.5 


3.4 x 10~ 3 




NO 


A 


0.5 


45 


8.8 


-1.72 


32.4 


3.1 x 10~ 3 


0.0821 


NO 






85 


10.7 


-0.75 


31.3 


3.6 x 10~ 3 




NO 






5 


3.9 


-1.99 


31.0 


9.5 x 10~ 3 




NO 


13 


0.75 


45 


3.5 


-1.59 


31.0 


8.5 x 10" 3 


0.112 


NO 






85 


4.1 


-0.67 


31.7 


9.9 x 10" 3 




NO 






5 


2.28 


-1.97 


31.1 


2.0 x 10" 2 




NO 


C 


0.875 


45 


2.07 


-1.54 


31.2 


1.8 x lO- 2 


0.145 


NO 






85 


2.43 


-0.56 


32.3 


2.1 x 10~ 2 




NO 






5 


1.9 


-1.68 


31.5 


3.5 x 10" 2 




NO 


D 


0.94 


45 


1.7 


-1.27 


31.8 


3.1 x 10~ 2 


0.179 


NO 






85 


1.86 


-0.44 


32.9 


3.4 x 10~ 2 




YES 






5 


1.86 


-1.37 


32.2 


6.0 x 10~ 2 




NO 


E 


0.97 


45 


1.68 


-1.01 


32.5 


5.4 x 10~ 2 


0.213 


NO 






85 


1.83 


-0.21 


33.8 


5.9 x 10" 2 




NO 






5 


1.4 


-1.48 


32.6 


0.11 




NO 


F 


0.98 


15 


1.23 


-1.13 


32.9 


0.10 


0.245 


NO 






85 


1.29 


-0.26 


34.3 


0.10 




NO 



Note. - 



Columns same as in Table [T] 
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Table 3. Summary of MHD and SED for models with Tj/T e =10. 



run 


a. 


i 

[deg] 


(M> 
[lO^Meyr" 1 ] 


a 


log 10 L x 
[ergs -1 ] 


'/ 


VTD 


consistent 
w/ obs.? 






5 


90.8 


-1.37 


30.1 


5.4 x 10" 4 




NO 


A 


0.5 


45 


117.1 


-0.2 


31.4 


6.7 x 10~ 4 


0.0821 


YES 






85 


369.0 


1.38 


33.5 


1.6 x 10~ 3 




NO 






5 


38.3 


-1.05 


30.0 


1.6 x 10~ 3 




NO 


B 


0.75 


45 


50.2 


0.04 


31.8 


2.0 x 10" 3 


0.112 


YES 






85 


190.6 


1.49 


34.6 


6.9 x 10" 3 




NO 






5 


19.5 


-1.15 


31.4 


5.1 x 10" 3 




NO 


C 


0.875 


45 


23.6 


-0.07 


32.0 


6.2 x 10" 3 


0.145 


YES 






85 


41.4 


1.19 


34.2 


1.7 x 10~ 2 




NO 






5 


13.7 


-0.93 


31.8 


1.1 x 10" 2 




NO 


D 


0.94 


45 


15.2 


-0.05 


32.3 


1.1 x 10~ 2 


0.179 


YES 






85 


31.2 


1.17 


34.4 


2.5 x 10~ 2 




NO 






5 


13.6 


-0.40 


32.5 


2.6 x 10~ 2 




YES 


E 


0.97 


45 


14.3 


0.2 


33.1 


2.8 x 10~ 2 


0.213 


NO 






85 


26.5 


1.19 


35.2 


5.1 x 10~ 2 




NO 






5 


9.07 


-0.6 


32.7 


5.3 x 10~ 2 




YES 


F 


0.98 


15 


9.16 


-0.08 


33.2 


5.2 x 10" 2 


0.245 


NO 






85 


15.5 


1.12 


35.4 


9.0 x 10~ 2 




NO 



Note. - 



Columns same as in Table [T] 
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Table 4. The size of the emitting region at 230 GHz. 









T-JT a 


= 4 


Ti/T c 


= 3 


Ti/T c 


= 40 


run 


a* 


i 


V\ 


a 2 


fl 


CJ-2 


01 


C2 






5 


44.7 


44.6 


38.4 


38.0 


35.9 


35.9 


A 


0.5 


45 


34.2 


28.4 


28.9 


25.9 


32.7 


34.6 






85 


23.4 


20.7 


23.3 


24.9 


37.6 


33.9 






5 


38.8 


38.7 


35.4 


35.3 


34.6 


34.6 


B 


0.75 


45 


28.6 


25.9 


26.7 


23.7 


34.2 


30.2 






85 


20.7 


49.4 


20.8 


20.4 


36.9 


32.3 






5 


39.3 


39.2 


35.2 


35.0 


34.9 


34.8 


C 


0.875 


45 


29.3 


26.7 


26.4 


23.4 


29.3 


28.0 






85 


20.6 


49.4 


20.4 


49.7 


30.6 


30.4 






5 


37.2 


37.4 


32.2 


32.4 


28.0 


28.0 


D 


0.94 


45 


27.4 


25.3 


24.0 


24.5 


24.7 


23.9 






85 


49.2 


48.4 


48.5 


47.2 


27.4 


26.4 






5 


37.0 


36.9 


34.2 


34.4 


26.5 


26.4 


E 


0.97 


45 


27.4 


25.7 


23.4 


24.6 


23.6 


23.4 






84 


49.6 


48.2 


48.4 


47.3 


25.8 


24.6 






5 


36.5 


36.4 


34.3 


34.2 


26.4 


26.4 


F 


0.98 


45 


27.2 


25.2 


23.4 


24.2 


22.8 


24.7 






85 


48.8 


48.4 


46.5 


47.6 


23.8 


22.6 



Note. — The columns from left to right are: run ID, dimensionless 
spin of the black hole, inclination angle i, size of the emitting region in 
terms of standard deviation in the major (tri) and minor ((J2) axis for 
Ti/Tc = 4 (col. 4 and 5), T;/T c = 3 (col. 6 and 7), and T;/T c = 40 (col. 
6 and 9) in units of f ias. For a Gaussian m odel, the VLBI data require 
FWHM= 37±{l fias jPoeleman et al.ll2008ri , or a = W^tl ^ as - 
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Density Magnetic field Temperature 




10 20 30 40 10 20 30 40 10 20 30 40 




Fig. 1. — Disk structure in our best-bet model with a* = 0.94 (model D) and with T-JT C = 3. 
The number density, the magnetic field strength, and the electron temperature are shown 
in the left, middle, and right panel respectively. The axis scale units are GM/c 2 .The figure 
presents a single time slice. 
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Fig. 2. — SEDs computed based on a single time slice t = 1680GM/c 3 (thick line) (see 
Figure [T]for the distributions of the physical variables corresponding to the same time) along 
with the tim e averaged spectrum (thin line) of our best-bet model. O bservatio nal points are 



taken from: iFalcke et all Il998l . lAn et all I2005L 



Marrone et al 



2006 at radio, Genzel et al. 



20031 at NIR (1.65, 2.16, and 3.76 /im) and Baganoff et al J 12003 at X-rays (2-8 keV). The 



upper limits in the NIR band ar e taken from 



Me 



Schodel et al.l 120071 (8 .6 urn) and iHornstein et al 



flarin g state are from I Genzel et al 



ia fc Falckdl200ll (30, 24.5 and 8.6 /an), 



20071 (2 lira). The po ints in the NIR at 



2001 



20031 (1.65, 2.16, and 3.76 /im) , and iDodds-Eden et al. 



2009 (3.8 /im). An example of X-ray flare (Lx = 1 x 10 35 ergs l ) is taken from lBaganoff et al. 



- 26 - 



Emission at 230 — 690 GHz Emission at 2—10 fj. m Emission at 2—8 keV 




Fig. 3. — Maps showing the point of origin for photons in our best-bet model. We show 
the logarithm of the sum of the photon weights in each zone, which is proportional to the 
number of photons seen at 230 — 690 GHz (left), 2 — 10/im (middle) and 2-8 keV (right) band. 
The axis scale units are GM/c 2 . The figure presents a single time slice at t — 1680GM/c 3 . 
Note that color bands differ in scales. 
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Fig. 4. — SEDs computed for T;/T e = 1, 3 and 10 in left, middle and right panels, respectively 
and i = 5 deg, 45 deg and 85 deg in bottom, middle and top panel, respectively. Models A, 
B, C, D, E and F, have different a* (see Tables HJ El and [3]). Each SED is a result of 
averaging over 200 individual SEDs taken from last 500GM/c 3 of the GRMHD runs (4 runs 
and 50 dumps for each spin). Observational points and upper limits as in Figure El Models 
consistent with the observations are marked with thick lines. 
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Fig. 5. — Images of the accretion flow at 230 GHz for models with SEDs that are consistent 
with observations of Sgr A*. The images has been averaged over time and over four separate 
realizations of each model. Intensities are given in units of ergs _1 pixelsize _2 Hz _1 sr _1 , where 
the pixel size is 0.82 /xas. The images show inner 40 GM/c 2 . Left top panel shows model 
D with T;/T c = 3 at i = 85deg. Left middle and bottom panels show high spin models E 
and F, respectively, for T-jT e = 10 and i = 5deg. The right panels show models A, B, and 
D for Tj/T c = 10 and i = 45deg in the upper, middle, and lower panel respectively. The 
white circle marks FWHM=37 fias of a symmetric Gaussian brightness profile centered at 
the image centroid. 



